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The  goal  of  this  work  was  to  develop  computational  procedures  for  the  calculation 
of  the  nonlinear  optical  properties  of  molecules  and  to  use  these  procedures  to  study  the 
molecular  nature  of  these  effects.  Of  particular  interest  have  been  studies  of  the  third- 
order  effects.  In  the  following  sections,  a  short  report  will  be  given  on  the  status  of  this 
work. 


Methods 


For  an  isolated  molecule  in  a  static,  dc-electric  field  the  energy  may  be  written  in 


terms  of  the  field  strength,  F,  as 

E(F)  =  E(0)  -  n,F,  -  ia,F,F^  -^P*F,F,F^  -  ;iT*F,F,F.F,  -  • 


(I) 


In  this  expression  ay  are  components  of  the  molecular  polarizability  and  and  yyi^]  are 
components  of  the  first  and  second  hyperpolarizability,  respectively.  Summation  over 
repeated  indices  is  assumed.  By  truncation  of  equation  1  after  the  y  term,  the  molecular 
properties  can  be  obtained  by  calculating  the  energy  at  several  values  of  the  applied  field. 
This  is  the  finite-field  method  as  has  been  implemented  in  MOP  AC  6.0  and  GAMESS 
using  the  procedures  presented  in  J  Comp.  Chem.  ^ 

Two  major  problems  with  the  finite-field  method  are  it  only  provides  static  field 
values  and  it  is  subject  to  several  numerical  difficulties.  Clearly  to  compare  with 
experimental  laser  studies,  the  proper  frequency  dependent  polarizability  and 
hyperpolarizabilities  are  required.  The  method  chosen  to  include  these  eflfects  was  the 
TDHF  procedure  of  Bartlett  and  Sekino.2  The  procedure  followed  is  very  well  described 
by  Kama  and  Dupuis.^  This  method  all  solves  the  numerical  difficulties  with  the  finite- 


field  procedure.  The  only  disadvantage  of  the  TDHF  approach  is  that  does  not  include 
electron  correlation. 

Oligomer  calculations 

a)  polyacetylene: 

The  initial  studies  focused  on  using  the  finite-field  procedures  implemented  m  the 
MOP  AC  semiempirical  program  to  treat  very  large  systems  such  as  oligomers  and 
polymers.  The  first  system  examined  was  a  series  of  polyene  oligomers,  H(C2H2)nH,  for 
n  up  to  20.'^  The  a  and  y  results  compare  very  well  with  a  series  of  ab  initio  calculations 
of  Hurst,  Dupuis,  and  dementi  .  5  This 


system  has  recently  been  reexamined  with  _  I/jj //  n _ 

the  TDHF  method  to  calculate  the  frequency  . ~ 

dependent  values  as  a  test  of  the  program  ^  «>  mmv/, 

and  to  explore  the  importance  of  these  ^ 

effects.  The  figure  at  the  right  shows  the  \  ,  ,  , 

timurn 

calculated  values  of  a(-<B;(i>). 

For  y,  frequency  dependence  means  there  are  several  sets  of  data  that  can  be 


obtained.  We  have  calculated 
three  different  y  curves:  ythg(~ 

a©;©,®,®),  yefish(-2©;ffl,©.0), 

and  yidir(-o;©,o,-®)  and  the 
figure  at  the  right  shows  a  typical 
behavior  of  these  quantities  for 
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H(C2H2)3H. 

In  this  work  we  also  studied  the  behavior  of  a  and  y  as  the  oligomer  length  is 
increased  by  adding  subunits.  Typical  behaviors  are  shown  below  for  a(-©;©)  and  y(- 
2ffl,©,©,0).  For  the  static  quantities,  a/subunit  seems  to  saturate  quite  quickly,  by  around 
12-14  units,  while  y/subunit  staturates  at  much  longer  oligomers,  about  20  units.  As  the 


frequency  increases,  the  values/subunit  still  seem  sto  staturate  but  at  a  higher  value  and  at 
longer  oligomers.  However,  if  the  frequency  is  increased  even  further,  the  quantities  start 
to  diverge  (i,e.,  no  saturation  behavior).  The  energy  at  which  this  seems  to  happen  is 


related  to  the  excitation  energy  of  the  polymer  (a  little  over  2  eV  for  polyacevylene). 


b)  Cyclopentadienes 

In  addition  to  the  polyene  oligomers,  the  heteroatom  substituted 
polycyclopentadiene  molecules  shown  below  were  studied. 

juuuyuUu  dMwyyuutf 
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These  systems  behave  similarly  to  the  polyene  oligomers  expect  they  saturate  the  per 
subunit  values  at  much  smaller  numbers  of  subunits.  These  systems  were  originally 
studied  using  the  finite-field  and  are  currently  being  reexamined  using  the  TDHF 
procedure. 

Hyperpolarizabiiity  Relationships 

General  formulas  for  third-order  nonlinear  optical  processes  have  been 

developed.^  One  such  useful  formula  is 

Y  =Y°(1  +  Ao)^+-)  , 


(2) 


where  .4  is  a  frequency  independent  constant  and  the  frequency  parameter  is  given  by 

The  importance  of  this  type  of  equation  is  that  if  it  can  fit  the 

data  when  truncated  at  the  A  term,  it  then  allows  the  calculation  of  the  dispersion  in  any 
type  of  Y  from  the  knowledge  of  a  fit  on  a  single  set  of  y  dispersion  data.  For  example  if  A 


is  determined  from  a  THG  experiment  (with  (t>l  =  12a>^),  then  equation  2  can  be  used  to 
calculate  the  dispersion  for  either  EFISH  (with  ml  =  6co^  or  DFWM  (with  ml  =  4©‘). 
When  valid  in  the  truncated  form,  this  equation  gives  the  simple  relationship  that  the 
dispersion  should  be  twice  as  great  for  THG  as  for  EFISH  and  three  times  the  amount  in 
DFWM.  If  this  type  of  expansion  is  valid  then  all  second  hyperpolarizabilities  would  be 
the  same.  Because  of  the  pole  behavior  of  f 

[  OKC 

y.  equation  2  can  not  be  totally  correct  but 

7  500 

it  does  work  well  for  small  frequencies,  as  ^ 

e  -too 

shown  to  the  right  where  the  results  for  ^ 

300 

CH2O  are  plotted  versus  ©l.  The  utility  of 

200 

relationships  between  different  types  of 
hyperpolarizabilities  demonstrates  an 
important  reason  for  adopting  the  convention  where  all  y  values  go  to  the  same  static 
limit. 

Vibrational  Corrections 

After  the  inclusion  of  frequency  dependence,  the  next  potentially  large  corrections 
to  the  calculations  of  hyperpolarizabilities  are  the  inclusion  of  vibrational  effects.  A 
dettdled  method  for  calculating  these  effects  has  been  worked  out  by  Bishop  and  Kirtman^ 
and  we  have  performed  detailed  SCF  and  MP2  calculations  on  water  as  a  test  molecule.® 
From  this  work,  the  estimated  vibrational  correction  for  P(0;0,0)  is  -37%  of  the  electronic 
value  and  for  y(0;0,0,0)  is  15%  For  the  hyperpolarizabilities  at  higher  frequency  (0.0656 
a.u.)  the  y  corrections  are  less  [y(-<o;©,0,0)  3%,  y(-2©;©,a,0)  -1%,  y(-3o;ffi;©;ffl)  -0%, 
and  y(-<a, ©,©,-©)  8%]  but  the  corrections  to  P  are  significant  [  3(-<a;©,0)  -19%,  P(- 
2© ;©.,©)  2%].  This  same  work  shows  however,  that  vibrational  corrections  have  the 
potential  to  be  enormous,  620%  for  P(0,0,0)  of  NH3. 


Intermolecular  Effects 

Another  effect  that  may  be  important,  particularly  in  condensed  phase  experiments, 
is  intermolecular  interactions.  Initially,  we  have  performed  ab  initio  calculations  on  the 
approach  of  two  ethylenes^  and  have  done  a  follow-up  study  for  two  acetylene  molecules 
The  main  results  from  these  studies  are  that  for  ab  initio  calculations  the  choice  of  basis 
set  is  critical  and  very  accurate  results  need  to  include  correlation  but  reasonable  results 
can  be  obtained  at  the  RHF  level. 

Recently,  we  have  been  studying  the  approach  of  two  benzene  molecules  The 
initial  calculations  have  been  done  at  the  AMI  level  and  show  some  interesting 
possibilities.  The  results  are  summarized  for  the  interaction  energy,  alpha,  and  gamma  for 
the  face-to-face  (circles),  side-on  (squares),  and  point-on  (diamonds)  approach  of  two 
molecules.  At  each  distance,  all  other  geometrical  parameters  where  optimized. 


The  above  results  show  very  little  effect  on  a  for  all  but  the  face-to-face  approach  while 
for  gamma  the  effects  are  much  larger,  especially  the  side-to-side  approach  where  the 


increase  in  y  is  about  150%.  The  energy  requirements  for  approaching  to  2.5  A  are  no 
more  than  about  40-50  kcal. 

Other  Studies 

Using  the  GAMESS  based  finite-field  procedure,  we  began  a  study  of  systems 
where  electron  correlation  effects  are  very  important.  Initial  calculations  were  performed 
on  the  Be  atom  because  the  near  degeneracy  between  the  and  Is22p2  levels  makes 

normal  perturbation  treatments  ill-behaved.  The  results  are  quite  good  and  are 
summarized  in  a  paper  in  Chem.  Phys.  Lett. 

References 

1.  H.  A.  Kurtz,  J.  J.  P.  Stewart,  and  K.  M.  Dieter,  J.  Comp.  Chem.  11,  82  (1990). 

2.  H.  Sekino  and  R.  J.  Bartlett,  J.  Chem.  Phys.  85, 976  (1986). 

3.  S.  P.  Kama  and  M.  Dupuis,  J.  Comp.  Chem.  12, 487  (1991). 

4.  H.  A.  Kurtz,  Intern.  J.  Quantum  Chem.  Sym.  24,  791-798  (1990), 

5.  G.  J.  B.  Hurst,  M.  Dupuis,  and  E.  Clementi,  J.  Chem.  Phys.  89,  385  (1988). 

6.  D.  M.  Bishop,  J.  Chem.  Phys.  90,  3192  (1989);  D.  M.  Bishop,  Chem.  Phys.  Lett  153, 

441  (1988);  D,  M.  Bishop,  7.  Chem.  Phys.  95,  5489(1991);  D.  P.  Shelton,  /  Chem. 
Phys.  84,  404  (1986);  Phys.  Rev.  A  37,  3813  (1988). 

7.  D.  M.  Bishop  and  B.  Kirtman,  J.  Chem.  Phys.  97 , 5255  (1992). 

8.  D.  M.  Bishop,  B.  Kirtman,  H,  A.  Kurtz,  and  J.  E.  Rice,  J.  Chem.  Phys.  98,  8024 

(1993). 

9.  D.  S.  Dudis,  A.  T.  Yeates,  and  H.  A.  Kurtz,  Mat.  Res.  Soc.  Symp.  Proc.  247, 93 

(1992). 

10.  T.  Pluta  and  H.  A.  Kurtz,  Chem.  Phys.  Lett.  189,  255  (1992). 


Personnel 


Below  is  a  list  of  all  the  personnel  at  Memphis  State  University  who  have 
participated  in  the  research  supported  by  this  grant. 


Dr.  Tadeusz  Pluta 

Dr.  WonH.  Kim 

Dr.  Farhad  Sayyarpour 
Pr^ashan  Korambath 
Shihyen  Chen 
Clayton  Spencer 

Shari  Jones 


Postdoctoral  Assistant  (currently  faculty  member  at  the 
Silesian  University,  Poland) 

Graduate  Student  (obtained  PhD  1991,  currently  faculty 
member  Cheju-do  National  University,  Korea) 

Graduate  Student  (obtained  PhD  1992) 

Graduate  Student  (graduating  Fall  1993) 

Graduate  Student  (MS,  1992) 

Undergraduate  Student  (from  Rhodes  College,  a  nearby 
liberal  arts  college) 

A  high  school  senior  (funded  by  a  university  sponsored  NIH 
program) 


